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Abstract. The reaction np — > ppn~ has been studied in a kinematically complete measurement with a 
large acceptance time-of-flight spectrometer for incident neutron energies between threshold and 570 
MeV. The proton-proton invariant mass distributions show a strong enhancement due to the pp( So) 
final state interaction. A large anisotropy was found in the pion angular distributions in contrast to the 
reaction pp — » pp^r - At small energies, a large forward/backward asymmetry has been observed. From 
the measured integrated cross section cr(np — > pp7r~), the isoscalar cross section aoi has been extracted. 
Its energy dependence indicates that mainly partial waves with Sp final states contribute. 
Note: Due to a coding error, the differential cross sections da/dM pp as shown in Fig. 9 are too small by a 
factor of two, and inn Table 3 the differential cross sections da/dfi* are too large by a factor of 10/2-7T. 
The integrated cross sections and all conclusions remain unchanged. A corresponding erratum has been 
submitted and accepted by European Physics Journal. 



1 Introduction 

Pion production is the basic inelastic process in nucleon- 
nucleon interaction. Triggered by the high precision data 
for proton-proton induced reactions [1] , a renewed interest 
arose within the last decade. It has been found that the 
existing theoretical description [2] at that time underesti- 
mated the near-threshold cross section data for pp — > d7r + 
by a factor of 1.8 [3] and for pp — * pp7r° even by a factor of 
5 [1] . Refinements by including the A isobar and the intro- 
duction of an energy dependence in the s-wave rescattering 
term provided an enhancement of the cross section predic- 
tion by about a factor of 2 [5] ■ Heavy- meson-exchanges [3] 
and the 'offshell' behaviour of the irN amplitude in the 
rescattering diagram [7] were discussed as possible mech- 
anisms and were both able to explain the discrepancy. 

In recent years, calculations with microscopic models for 
the NN and 7rN interactions based on meson exchange 
were developed by groups in Jufich 8.9, lOjTT] and Os- 
aka [T31|T3] . At present, these are the only models consid- 
ering all single pion production channels including higher 
partial waves. 

In addition, first calculations in chiral perturbation the- 
ory were performed for neutral pion production, 14, 15, 16, 
I17j and charged pion production [18j . For neutral pion pro- 
duction also one-loop diagrams in the formalism of heavy 
baryon chiral perturbation theory were calculated [19,20, 



[21] . Recently, also higher partial waves were calculated 
in the framework of chiral perturbation theory |22j . It was 
also shown, that the information that can be deduced from 
pion production in nucleon-nucleon collisions is relevant 
for constraining three nuclcon forces. In the context of 
this presentation it is important that the relevant oper- 
ator structure that allows this connection contributes to 
the reaction np — > pp7r _ as well. 

To pin down the different production mechanisms, high 
precision data from different pion production reactions are 
needed. This paper addresses an improvement of the data 
quality for charged pion production in neutron-proton col- 
lisions in order to achieve a better knowledge on the isoscalar 
cross section ctqi- 
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2 Pion production in np collisions 

Under the assumption of isospin invariance, all single pion 
production reactions in nucleon-nucleon collisions into a 
three body final state can be decomposed into three par- 
tial cross sections a^i, , where I{ and If denote the isospin 
of the two-nucleon system in the initial and final state, re- 
spectively [23] (see Tab.[T]). At medium energies, the cross 
sections an and a w are dominated by the excitation of 
the intermediate ^33 resonance and are well-measured 
even close to threshold [1,24.25,26,27 . In contrast, the 
isoscalar cross section aoi, which has to be extracted from 
pion production data in neutron-proton and proton-proton 
collisions (see Sec. l2.lj) . is not well-known. Due to isospin 
conservation, the NZ\ intermediate state is not accessible 
from an I = initial state and therefore eroi is expected 
to be small if resonance production dominates. 



2.1 Determination of <r i 

The cross section aoi can be extracted from np — > NN7r ± 
data by measuring 
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is obtained. The an cross section in the intermediate en- 
ergy range is well-known from pp — > pp7r° measurements. 
However, the situation for cr i was not clear in the past. 
Several experiments reported significant coi values be- 
low 600 MeVpgpgilBUllBTP^. while others found small 
or even negligible aoi contributions for energies up to 750 
MeV [33l[341135"] . A partial wave analysis of Arndt and Ver- 
west 36[ gave no significant er i contribution below 1 GeV 
while Bystricky et al. [37 found small, but non-negligible 
values in a similar analysis. This unsatisfactory finding 
may be addressed to large experimental uncertainties and 
inconsistencies in both, the pp and the np data at that 
time. The determination of ctoi from several former np 
cross section measurements suffered from averaging the 
results over a large neutron beam energy range [58 , 30 , 
[38] . The pp — * pp7T° data have been remarkably improved 
during the last decade. Hence, the energy dependence of 
°np^NN7r± has to be determined with much higher preci- 
sion than it was obtained by former experiments in order 
to extract troi reliably. 

The determination of am includes a principal model de- 
pendence as © only holds in the case of exact isospin 
invariance. However, due to the different particle masses 
entering into the reactions np — > pp7r~, np — > mi7r + and 
pp — > pp7T°, isospin invariance is only an approximate 
symmetry. As a consequence, the comparison of the cross 
sections can not be performed at the same beam energy. 
Two methods have been discussed in the literature [3TJ 
I3"2"] so far which will be denoted in the present discus- 
sion as the r\- and the y^-scheme, respectively. In the 
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Table 1. Decomposition of pion production cross sections in 
partial cross sections cr^i, . 



77-scheme, the subtraction is performed at equal values of 
7] = p* max /m w + , the maximum value of the dimensionless 
cm. pion momentum. In the y^-scheme the two reactions 
are compared at the same cm. energy yfs. This corre- 
sponds to a resonant production mechanism with a AN 
intermediate state. However, near the production thresh- 
old, it neglects the different threshold values for the three 
reactions. 

A modification of the -^/s-scheme, the Q-scheme, performs 
the subtraction at equal excess energies Q = ^/s — y / sth7 
above the cm. threshold value y / Sthr- In the -y/s-scheme 
the differences in beam energy are AT = T n — T p = 
—2.6 MeV, whereas for the 77- and the Q-scheme, the dif- 
ference in beam energy is quite similar and reads AT w 
+7 MeV at T n = 350 MeV and AT « +4 MeV at T n = 
550 MeV. Since the cross section rises strongly between 
threshold and about 700 MeV beam energy, the results on 
<7oi in the y^s-scheme on one hand and in the 77- and the 
Q-scheme on the other hand differ significantly. 



2.2 Angular distributions 

A different approach to establish the existence of <toi takes 
advantage of the properties of the pion angular distribu- 
tions. All single pion production amplitudes with three- 
body final states can be decomposed in terms of three 
isospin amplitudes Mj f j i [37] which are related to the par- 
tial cross sections aij f by 
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The amplitudes for the pion production reactions of inter- 
est then become 

<pp7r°|A/jpp> =-<nn7r°|M|nn>= -^A'/n (6) 

V2 



<pp7r |M|pn > = — < mi7r |M|np > = 
<pp7r~|M|np> = -<mi7r + jM|pn> = 



J_Mi + iMu (7) 
J_Mio + iMn. (8) 
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Some consequences follow from these relations. Due to 
the identical particles in the initial state, the pion cm. 
angular distribution in the reaction pp — * pp7r° is for- 
ward/backward (f/b)-symmetric. Hence, a f/b-asymmetry 
observed in the reaction np — » pp7r~ or np — » nn7r + in- 
dicates the presence of cr i caused by an interference be- 
tween the amplitudes M\q and M\\. The same conclusion 
holds if differences in the cross sections at the same pion 
cm. angle are found for the reactions np — ► pp7r~ and 
np — > nn7r + . 

The differential cross section for the reaction np — * NN7r ± 
can be expanded in terms of powers of cos 9* |39j : 

^ = ao ± ai • cos 0* + a 2 • cos 2 9* ± (9) 

where the '+' and '— ' sign corresponds to the charge of 
the pion. Assuming that pion orbital angular momenta 
£ > 1 can be neglected, the expansion is truncated af- 
ter the quadratic term. A f/b-asymmetry in the reaction 
np — > NN-7T is then described by a non-vanishing linear 
coefficient oi. Most experiments measuring the reactions 
np — > NN-zr^ found significant linear cosine terms for en- 
ergies well below 600 MeV [28lf39~l and small or vanish- 
ing values above 600 MeV [3~3l|381130"] . No significant f/b- 
asymmetry was reported by Bachman et al. [40] at 443 
MeV. 

Historically, pion angular distributions in proton-proton 
reactions were parametrized by 

+ ,io, 

where b is called the anisotropy parameter. To include an 
angular asymmetry for the reaction np — ► NN-zr^ , one may 
add a cos 9* term resulting in 

^ =C- Q + a-cos^ + fo-cos 2 ^ . (11) 

This parametrization is still appropriate for comparison 
with older data. Below 600 MeV, most experiments found 
anisotropy parameters b for the reaction np — > NNzr* 
which were significantly larger than those of the reaction 
pp — > pp7T° [28,31,39,40] indicating the presence of the 
(Tax cross section in np — > NNzr*. Nevertheless, no con- 
clusion concerning the energy dependence of 6 n p^NN7r± 
can be drawn from the existing data sets (see Sect. l5.2|) . 
since they are not fully compatible. Moreover, some ex- 
periments [31 , 39 were restricted by acceptance cuts and 
could extract the parameters only in a model-dependent 
way. 

2.3 Partial waves 

A helpful tool for the understanding of the production 
mechanism is provided by a partial wave decomposition 
of the scattering matrix. In a usual coupling scheme, the 
partial wave is written as 23+1 Lj —t 2S +1 L'j,£j, where S is 



the total spin, L the orbital angular momentum and J the 
total angular momentum of the two nucleons in the intial 
state, while S' , L' and J' give the corresponding angular 
momenta in the final state. The orbital angular momen- 
tum of the pion with respect to the final state di-nucleon 
system is denoted by I. The conservation of angular mo- 
mentum, isospin and parity and the consideration of the 
Pauli principle for the di-nucleon system in the inital and 
final state lead to a remarkable reduction of possible par- 
tial waves. In particular, in the reactions np — > NN-zr^ 
and pp — > pp7r°, partial waves with Ss final states can 
only contribute to u\\ and partial waves with Sp final 
states only to a m . There is only one possible Ss partial 
wave, 3 Po — ^SoSo, whereas there are two possibilites that 
lead to Sp partial waves, 3 Di — ^SoPi and 3 Si — *- 1 SoPi- 
Each partial wave shows a characteristic angular depen- 
dence as a function of 9*, which is constant in the case 
of 3 Po — > 1 SoSo and 3 Si — ^SoPi while it is described by 
| + cos 2 9* for 3 Di — > 1 SoPi- A further interesting feature 
of partial waves is the expected 77 dependence of their ex- 
citation function. If the pp final state interaction can be 
neglected, as in the case for pp P-waves, an 77 dependence 
oc rj 2 '( L + £ + 2 ) is expected [41]. In contrast, for pp S-waves 
the final state interaction plays an important role. In this 
case, an excitation function of the form oc is ex- 

pected [41]. However, in the reaction pp — > pp7r° close to 
threshold, where only the partial wave 3 Po — ^SoSo con- 
tributes, a clear deviation from the naive ^-dependence 
was observed [JJ . 



2.4 Recent developments 

Considerable improvement has been achieved during the 
last decade with new medium-energy-accelerators which 
provided secondary neutron beams of high intensity and 
high polarisation as well. Below the two-pion production 
threshold, single spin observables in the reaction np — * 
pp?!- - have been measured at TRIUMF at 443 MeV [40] 
and at SATURNE at 572 MeV [35]. Exclusive experiments 
at TRIUMF with proton beams at energies of 353, 403 and 
440 MeV incident on a deuterium target were dedicated 
to investigate the pp( 1 So) final state [43;,,44j. They revealed 
the significance of the itqi cross section in that particular 
phase space configuration by a subsequent partial wave 
analysis considering the partial waves 3 Si ^ x S0P1 and 
3 Di -^SoPi for the I = and 3 P -^SqSq for the 1 = 1 
initial state [33J • At 440 MeV, even a small contribution 
from pion d- waves 3 P2 — > 1 Sod2 and 3 F 2 — s- 1 Sod2 was 
reported. 

Despite of this progress, the data of pion production in 
neutron-proton collisions are still lacking precise measure- 
ments of differential and integrated cross sections. In this 
paper, we report on a kinematically complete measure- 
ment of the reaction np — * pp7r~ using a polarised neutron 
beam. Here, we deal only with spin averaged results while 
spin dependent observables will be presented in a forth- 
coming publication. Details of the analysis can be found 
ingS]. 
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3 Experiment 

3.1 Neutron Beam 

We used the polarised neutron beam facility NA2 at 
Paul-Scherrer-Institut (PSI) which is described in detail 
in Ref. [46]. Vertically polarised protons from an atomic 
beam source were accelerated in the cyclotron to an en- 
ergy of 590 MeV. The beam polarisation vector was then 
rotated from vertical into longitudinal direction and was 
reversed every second at the source. Longitudinally po- 
larised neutrons were produced in the reaction 12 C(p, n)X 
on a 12 cm thick carbon target [4"6"ll47] . Neutrons emitted 
at 0° with respect to the proton beam axis were selected 
by means of a collimator of 2 m length. The neutron beam 
was stabilized on its axis using a feed back system which 
kept the proton beam at the center of the neutron pro- 
duction target within 0.1 mm [48j . Remaining beam pro- 
tons and secondary charged particles were deflected by a 
dipole magnet. A lead filter reduced the 7-contamination 
in the beam, originating mainly from the decay of neutral 
pions which were produced in the neutron production tar- 
get. Spin rotating magnets allowed to choose any neutron 
polarisation direction. All results presented in this paper 
were achieved by averaging the data for the two beam po- 
larisation states for a transversly polarised neutron beam. 

The time structure of the PSI proton beam consists of 
bunches of 0.84 ns width (FWHM) with a bunch frequency 
of 50.63 MHz. The neutrons show a continuous energy dis- 
tribution with a quasi-elastic peak at about 530 MeV and 
a broad continuum at lower energies resulting mainly from 
pion production and A excitation [49] . For a typical proton 
beam current of 10 /iA and a beam collimator opening of 
9 mm diameter, a neutron beam with a typical flux of sev- 
eral 10 7 n/s was obtained. The beam was 4 cm in diameter 
(FWHM) at 20 m downstream of the neutron production 
target. 

Two monitors, described in [46| , were used to check the 
neutron beam properties during data taking and in the 
offline analysis. Monitor I was placed immediately behind 
the neutron beam pipe window in the NA2 area and al- 
lowed a relative measurement of the neutron beam in- 
tensity using the reactions H 1 (n, p)n and C(n,p)X in a 
Polyethylene target. It consisted of three scintillator coun- 
ters Ml, M2 and M3 with the Polyethylene target sand- 
wiched between Ml and M2. The counting rate Ml- Ml ■ 
M3 served as a measure of the neutron beam intensity. 
Monitor II was placed 2 m behind the experimental set- 
up. It measured the beam position and the beam profile 
using a scintillator hodoscope. The beam polarisation in 
horizontal and transverse direction was monitored by two- 
arm polarimeters. For the two different polarisation states 
at the source (flipped and non-flipped) , the neutron beam 
intensities, positions and polarisations were found to be 
equal within the statistical errors. 

3.2 Experimental set-up 

For the kinematically complete measurement of the re- 
action np — > pp7r~, a time-of- flight (TOF) spectrometer 



with a large geometrical acceptance was used. In the re- 
action np — > pp7r~, both protons are emitted within the 
angular range 0° < 8 p < 45° in the laboratory system for 
neutron kinetic energies below 570 MeV. The pion emis- 
sion angle is not constrained for any neutron energy above 
313 MeV. The experiment relied on the determination of 
the energy for each incident neutron provided by a TOF 
measurement and the reconstruction of the emission an- 
gles and velocities of at least two of the three charged 
particles in the final state. The experimental set-up (Fig. 
[IJ consisted of a liquid hydrogen target, two sets of drift 
chambers (DC8, BD6) as tracking devices, a segmented 
trigger hodoscope and a large area TOF wall. 



TOF wall 

BD6 




m 1 m 2 m 3 m 

Fig. 1. Top view of the detector. 



The lens-shaped target cell with walls consisting of a 125 
/im thick Kapton layer was filled with liquid hydrogen and 
placed 20.03 m downstream of the neutron production tar- 
get. It was 9.3 cm in diameter and 3.15 cm thick at central 
incidence and surrounded by several layers of superinsu- 
lation with a total thickness of 340 /im. The target cell 
was placed inside a vacuum vessel with an entrance and 
exit window of 10 cm and 30.8 cm diameter, respectively. 
Both were sealed by titanium sheets of 25 fim (entrance) 
and 50 (im (exit) thickness. 

The drift chamber DC8 with an active area of 56 x 56 cm 2 
was placed immediately behind the target vessel. It con- 
sisted of eight planes with alternating wire orientations 
in vertical (Y) and horizontal (X) direction. Each plane 
contained 14 cells of 4 cm width equipped with pairs of 
signal wires with a wire spacing of 0.42 mm in order to 
avoid left-right ambiguities for a single track crossing the 
drift cell. The large drift chamber stack BD6 consisted of 6 
planes with an active area of 214 cm(X) x 114 cm(Y) and 
with wire orientations UYVUYV where U and V denote 
directions of ±30° with respect to the Y direction. Each 
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drift cell had a width of 2 cm and was equipped with one 
signal wire. 

For both drift chambers, a gas mixture of 67.8 % Argon, 
29.5 % Isobutane and 2.7 % Methylal was used. The aver- 
age detection efficiency per plane extracted from one-track 
events (mainly protons) was found to be 97% for DC8 and 
95% for BD6. The spatial resolution for the DC8 planes 
was determined to be a = 0.2 mm for track angles of 0° 
with respect to the neutron beam axis increasing up to 
0.5 mm at 40° in agreement with former findings for the 
same chamber type [5D] . 

The trigger hodoscope which was placed immediately 
behind the DC8 had an active area of 65 x 65 cm 2 and 
consisted of two planes with horizontal(h) and vertical(v) 
scintillator slabs. Both planes were built up by twelve 
scintillator slabs with an area of 65 x 5 cm 2 and a thick- 
ness of 0.3 cm. Each slab was instrumented from both 
sides with HAMAMATSU R1450 photomultiplier tubes. 
In the central part of both planes, two scintillators (30 
x 5 x 0.3 cm 3 ), instrumented with only one photomulti- 
plier tube, were placed in order to leave a quadratic gap 
of 5 x 5 cm 2 for the neutron beam. By this gap, genera- 
tion of a large background rate in the detector by beam 
interactions inside the hodoscope was avoided. For a fast 
timing BICRON BC404 scintillator material was chosen 
in conjunction with BC800 as lightguide material in order 
to minimize light losses in the UV range. 
Three meters downstream from the hydrogen target cell a 
large scintillator TOF wall with an area of 300 x 290 cm 2 
was installed. The wall consisted of two groups of seven 
BC412 scintillator bars with dimensions of 300 x 20 x 
8 cm 3 instrumented on both sides with Philips XP2040 
photomultiplier tubes. The bars and photomultipliers had 
been already used in the LEAR experiment PS199 [51j. 
The two groups were separated by two smaller BC408 
scintillator bars with dimensions of 100 x 10 x 10 cm 3 
viewed from the outer side by Philips XP2020 photomul- 
tiplier tubes. They were placed to leave a gap of 10 x 10 
cm 2 in the centre of the TOF wall for the neutron beam. 



3.3 Electronics and data acquisition 

The electronics for all scintillator detectors had been 
already successfully operated in the LEAR experiment 
PS199 [51] and the PSI experiments on elastic np scatter- 
ing [52]. Compact CAMAC-modules, constructed by the 
University of Geneva, with constant-fraction discrimina- 
tors and mean-timers (DPNC 982) provided output sig- 
nals for ADC and TDC measurements as well as for trig- 
ger building. The TDC measurement was performed with 
timc-to-charge converter (TQC) modules, constructed by 
the University of Geneva, and LeCroy 4300 ADC CAMAC- 
modules for charge digitization, resulting in a time reso- 
lution of 50 ps. For the ADC measurement, LeCroy 4300 
CAMAC-modules were used. A dedicated ADC and TDC 
channel readout by the data acquisition software accord- 
ing to the fired scintillators was realized using LeCroy 
4448 coincidence register modules. The number of hits 
for each scintillator plane was determined by a LeCroy 



4532 majority logic unit (MALU) generating output sig- 
nals in case of at least one hit (ORO) and at least two 
hits (MDO). A traversing of the hodoscope by at least 
two charged particles was indicated by the coincidence 
signal TCP = (ORO h A MDO v ) V (MDO h A ORO v ). The 
subscripts h, v refer to horizontal and vertical orientation, 
respectively. 

Both drift chambers were instrumented with electronics 
cards, developped by the University of Freiburg, perform- 
ing pre- and main-amplification as well as discrimination 
of the drift chamber signals [53] . Drift time measurement 
was provided by LeCroy 1879 FASTBUS modules with a 
TDC resolution of 2 ns. Moreover, the hit multiplicity for 
each DC8 plane was determined using multiplicity cards, 
constructed at the University of Freiburg and successfully 
operated at the LEAR experiment PS202 [54]. The output 
signals of these cards were used in the second level trigger 
stage. 



hit 
multi- 
plicity 











TOF __ 

wall 

DAQ ready ■ 
50 MHz ■ 



DPNC 

982 



MALT" 



= fT" A 

► > reset 



TTOnOSCOPF. 



DPNC 
982 



DPNC 
982 



MALT 



2 nd level 
rigger 



com men 
» start /stop 
gate ADC 



► |l:200| -» 



spin 
flip 



C spin flip 
event 



scaler 
read out 



Fig. 2. Overview of the trigger electronics. 



Three different types of triggers were used in the exper- 
iment. A simplified scheme of the trigger electronics is 
shown in Fig. [2] 

A) The main trigger consisted of two levels and intended 
to select events of the reaction np — + PPtt~ by requir- 
ing a charged multiplicity of at least two. The first level 
trigger was built by a coincidence of i) the 50.63 MHz ra- 
dio frequency signal of the accelerator (50 MHz), ii) the 
TCP coincidence signal and iii) a computer- ready signal 
from the data acquisition system. The 50 MHz signal indi- 
cated the arrival time of the proton bunch at the neutron 
production target modulo 19.75 ns which is the cyclotron 
repetition time. It determined the timing of the first level 
trigger which served as a common stop signal for the drift 
chamber FASTBUS TDCs as well as a common start sig- 
nal for the hodoscope and TOF wall scintillator TQCs 
providing the times TOFhodo and TOF wa x\. Both times 
are measured modulo 19.75 ns. They include the time- 
of-flight from the neutron production target to the liquid 
hydrogen target (TOF\) and in addition the short time- 
of- flight to the hodoscope (TOF2) and respectively to the 
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TOF-wall {TOF 2 + TOF z ), (see Fig.©: 

TOF hodo = TOF 1 + TOF 2 + i ■ 19.75 ns (12) 
lW W aii = TOF 1 + TOF 2 + TOF 3 + i ■ 19.75 ns (13) 

where i = 0, 1, .... 



category 


measured observables 


CI 


2-prong + 2 velocities 


C2 


3-prong + 1 velocity 


C3 


3-prong + 2 velocities 


C4 


3-prong + 3 velocities 



Table 2. Fit categories and measured observables. 
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Fig. 3. Principle of the time-of- flight measurements. 

The second level trigger required at least one hit in 
the TOF wall (Oi?O wa n) and at least one hit in each 
DC8 plane. The latter condition suppressed a large part 
of events originating from beam interactions in the ho- 
doscope or the DC8 chamber. If the second level trigger 
condition was not fulfilled, a fast-clear signal was sent to 
all modules. 

B) A minimum bias trigger which required, at the first 
level, a minimum charged multiplicity of one in the ho- 
doscope and no further requirements at the second level 
was accepted with a prescaling factor of 1:200 correspond- 
ing to about 20% of the total trigger rate. These events 
were mainly elastic np scattering events giving by far the 
largest contribution to the total np cross section in this 
energy range. 

C) A pseudo-event trigger signal was generated at the 
end of each one-second-period and all main and minimum 
bias trigger signals were rejected during a 10 ms INHIBIT 
signal. During this time interval the proton beam polarisa- 
tion at the beam source was reversed and scaler registers 
were read out and reset. 

The data acquisition and readout was controlled by a 
STARBURST ACC2180 front end computer which buffered 
events and sent them in data packages to a VAX4090 
workstation where they were written to tape. Typical trig- 
gered event rates were of the order of 260 events per sec- 
ond with an average event length of 250 Bytes. The life 
time of the data acquisition system was about 30 — 40%. 
For the determination of the spin averaged results, about 
6 • 10 7 events with transversely polarised neutrons were 
recorded. This includes also dedicated calibration runs. 
About 2 • 10 7 triggers with transversely and longitudinally 
polarised neutrons were taken with an empty target cell in 
order to study background contributions originating from 
the target surroundings. 



4 Event reconstruction 

The reconstruction of the reaction np — > pp7r~ relied on 
the energy determination of each incident neutron by a 



time-of-flight measurement and the measurement of a suf- 
ficient set of track parameters and velocities for the emit- 
ted particles. For a given neutron energy, nine kinemati- 
cal observables, e.g., the momentum vectors of the three 
particles in the final state, describe the kinematical state 
completely. Due to energy-momentum conservation, it is 
sufficient to measure five of those. Events of the reaction 
np — > pp7r~ were reconstructed using a kinematical fit of 
at least six measured kinematical observables, and thereby 
separated from background. As a consequence, the reac- 
tion np — > pp7r _ could only be reconstructed from 2- and 
3-prong events but not from 1-prong events. However, 1- 
prong events from the minimum bias trigger sample were 
used to reconstruct elastic np scattering events. In the 
case of 2-prong events, the velocity for both tracks had 
to be measured in order to perform a kinematical fit. For 
3-prong events, the six track angles, in principle, were al- 
ready sufficient. Nevertheless, at least one particle velocity 
was required in order to determine the neutron time-of- 
flight (see Sec. 14.3ft . The particle velocities were given by 
two different types of time measurements. In general, they 
were determined from the time-of-flight TOF3 between the 
hodoscope and the TOF wall. In some cases, they were 
given by a determination of the time-of-flight TOF2 be- 
tween the liquid hydrogen target and the hodoscope (see 
SecSU). 

If one of the measured tracks could be associated to a 
pion or two measured tracks could be assigned to protons, 
the final state was completely fixed. This association was 
performed using kinematical arguments. Otherwise, up to 
three configurations were possible. In this case, the con- 
figuration with the best x 2 defined the particle-to-track 
association. 

According to the number of measured velocities in 2- 
and 3-prong events, four categories were defined for the 
kinematical fit. Their classification, CV, corresponding to 
the number of degrees of freedom v, is shown in Tab.O 



4.1 Track reconstruction 

Tracks were searched using the drift chambers BD6 and 
DC8. In the first stage, three-dimensional track intersec- 
tion points were reconstructed inside the large drift cham- 
ber (BD6) by combining only the positions of the hit wires 
in the U-, V- and Y-planes. This procedure speeds up the 
track finding and gives the BD6 intersection point with a 
sufficient precision of approximately ±20mm for the sub- 
sequent track finding in the DC8. 

In the second stage, track intersection points in the BD6 
and hits in the first DC8 X- and Y-plane provided esti- 
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mates for the track parameters. Hits were then searched in 
the other DC8 planes inside defined corridors. If at least 
three hits in both, the X and Y projection, could be found, 
the track candidate was accepted and a straight line was 
fitted to these DC8 hits. If more than one candidate was 
found for a given BD6 intersection point, only the can- 
didate with the best confidence level in the track fit was 
kept. Track candidates which gave no intercept with the 
liquid hydrogen target cell were rejected from the further 
analysis. 
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4.2 Scintillator Information 

Scintillators were associated to tracks if i) a track in- 
tercept with the fired scintillator was found and ii) the 
scintillator coordinate, calculated from the time difference 
between both photomultiplier signals using the measured 
(effective) scintillator light velocity, corresponded to the 
track coordinate at the scintillator. Light velocities and 
spatial resolutions were found to be 0.47 • c and 20 mm 
(cr) for the hodoscope and 0.53 • c and 45 mm (cr) for the 
TOF wall corresponding to time resolutions of 0.135 ns 
for the hodoscope scintillators and 0.280 ns for the bars. 

Time calibrations for scintillators were performed in sev- 
eral steps. Special runs without lead filters provided a 
large fraction of high energetic photons in the beam orig- 
inating from decays of neutral pions in the neutron pro- 
duction target. A part of these photons converted in the 
LH2 target and its surroundings into relativistic electron- 
positron pairs which gave a unique time signal for both 
times, TOFhodo and TOF wa n. Since the conversion process 
favours strongly the forward direction, only hodoscope and 
TOF wall scintillators near the beam line could be cali- 
brated reliably by this method. All bars in the TOF wall 
far from the beamline were calibrated with events from 
the minimum bias trigger sample requiring the elastic np 
scattering kinematics to be fulfilled. The remaining ho- 
doscope scintillators for a given plane were calibrated by 
means of the well-calibrated inner scintillators of the other 
plane using tracks that had traversed both scintillators. 



Fig. 4. Histogram: z-positions of reconstructed vertices z ve rtex 
for 2- and 3-prong events. For the other distributions see 
Sec. \4-6\ Hatched: events with a confidence level CL > 0.05 
in the kinematical fit. Cross-hatched: events of 'pit' -type with 
CL > 0.05. Horizontal lines: events of 'pn'-type with CL > 
0.05 for data with an empty target cell. 



Finally, a corridor track search inside the DC8 was per- 
formed using only those hits which had not been associ- 
ated to any track before. For this purpose, it was required 
that any new track candidate should intersect the liquid 
hydrogen target cell together with other tracks which had 
been found already in the second stage. This allowed to 
find tracks for particles which indeed traversed the DC8 
but stopped or, in case of a pion, decayed eventually before 
reaching the BD6. 

For 2-prong events, the vertex was defined as the point 
of closest approach while for 3-prong events, the vertex 
of the event was defined as the average of the three 2- 
track vertices. A typical vertex resolution of the order of 
3 mm was found. The z-position of the reconstructed ver- 
tices in the target region (Fig. [4]) shows a clear distinc- 
tion between events originating from the target cell and 
its surroundings on one hand and events from the tar- 
get entrance and exit window and the DC8 entrance win- 
dow on the other hand. For the further analysis a part of 
background events was rejected by a loose vertex cut of 
70 mm < Zvcrtcx < 135 mm. 



4.3 Neutron time-of-flight 

The time-of-flight TOF n for a neutron with velocity (3 n 
was given by the time TOFhodo — TOF2 where TOF2 was 
determined from the particle velocity Pi (see Sec. l4.4[) . For 
2- and 3-prong events, the track with the fastest velocity /3i 
was used for the determination of TOF n in order to min- 
imize systematic errors due to the energy loss correction 
(see Sec. l4.4p . In general, the time TOF n was ambiguous 
since the time TOFhodo was measured only modulo 19.75 
ns. However, due to the reaction threshold of 287 MeV 
neutron kinetic energy and the chosen distance of 20.03 
m between the production target and the liquid hydro- 
gen target, TOFhodo was unambiguous for np — > pp7r~ 
events. Nevertheless, two values for TOF n were possible 
since TOF2 was slightly different for pions and protons 
(see Sec. l4.4|) . The time resolution of TOF n was dominated 
by the width of the incoming proton bunch of 0.355 ns (cr) 
resulting in an uncertainty of 3 MeV (cr) at 287 MeV neu- 
tron energy and 9 MeV (cr) at 570 MeV. 

4.4 Velocity determination 

The velocity /3j ncas along the flight path between the ho- 
doscope and TOF wall scintillator was determined from 
the time-of-flight TOF 3 = TOF wall -TOF hodo . From /3f cas 
the initial velocity (3j was calculated with an energy loss 
correction function depending on the track angle and the 
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track origin inside the liquid hydrogen target using the 
GEANT 3.21 package [56]. For 2- and 3-prong events, the 
track origin was defined by the vertex while for 1-prong 
events it was assumed to be the track intercept with the 
centre plane of the target cell. The correction also depends 
on the particle type, and hence two possible velocities, /3f 
and f3f , have been associated to the track. Protons with 
(3f < 0.28 and pions with /3f < 0.45 stopped before reach- 
ing the TOF wall. The relative uncertainty 8j3i/ j3\, domi- 
nated by the time resolution of the TOF wall scintillators, 
was 3.3% at A = 1 and 1% at ft = 0.3. 

For several reasons, it was not always possible to de- 
termine the velocity (3\ for all tracks in a given event: i) 
several tracks had traversed the same scintillators, ii) the 
emission angle of a track could lie outside the TOF wall 
acceptance, iii) a particle stopped or iv) a pion decayed 
before reaching the TOF wall. This lowered the recon- 
struction efficiency for 2-prong events substantially. How- 
ever, the problem was circumvented if (3\ was measured 
for one track, denoted by track 1, and at least TOFhodo 
was available for the other track (track 2). In this case, the 
timc-of-flight for track 2, 

TOFi 1 = TOFfodo (track 2) - TOF n (track 1), (14) 

was determined resulting in a velocity measurement, de- 
noted /3JY cas - The error due to the time uncertainty of the 
incoming proton bunch cancels in (fl4|) . As a consequence, 
the uncertainty on TOF^ 1 was given by the timing resolu- 
tion of the hodoscope scintillators traversed by track 2. To 
reduce the uncertainty on /3JY° as , only tracks with /3JY° as in- 
formation in both hodoscope planes were accepted. This 
resulted in the relative uncertainty 5/3n//?n of 8.2% for 
relativistic particles. 




Fig. 5. Deposited energy in a TOF wall scintillator bar as a 
function of 7bar- Dashed line: expectation for stopped pions. 
Signals in the white region were considered as originating from 
protons for the kinematical fit. 



only about 1:6 from the reaction np — > pp7r~ is expected 
within the TOF wall acceptance where most pions have 
large momenta; ii) the probability for a pion decaying be- 
fore reaching the TOF wall is more than 50% for 7 < 1.25; 
iii) pions with 7 < 1.12 stopped before reaching the TOF 
wall. For the particle-to-track association, all tracks with 
entries outside the grey area could not be produced by 
a pion and hence were considered as originating from a 
proton. 



4.5 Particle-to-track association 

No dedicated particle identification (PID) was forseen in 
the experiment since the detector was designed as a TOF 
spectrometer. Therefore, the PID was obtained from the 
kinematical fit procedure. However, if for a given neutron 
energy the particle velocity superceeded the maximum 
possible value for a proton in the reaction np — ► pp7r~ 
by three standard deviations, the track was assumed to 
belong to a pion. 

In addition, the combined information from the deposited 
energy -Edep in the TOF wall and the 7-factor 7b ar = 
l/^/l — 0\ bar at the scintillator bar was used in order 

to exclude the pion hypothesis for certain tracks (Fig. 
[5]). The deposited energy -Edep increased with decreas- 
ing 7bar due to the energy loss described by the Bethe- 
Bloch formula. Pions with 7t, ar < 1-33 and protons with 
7bar < 1.1 stopped and deposited all their kinetic energy 
Ikin = m ■ (7bar — 1) inside the TOF wall. Clear signals 
from stopped protons and deuterons are seen. The ob- 
served deviation from linearity in the data is caused by 
a known saturation effect of the meantimer-discriminator 
module [5T] . No significant signal from stopped pions is ob- 
served due to several reasons: i) a pion-to-proton ratio of 



4.6 Kinematical fit and background rejection 

The events were reconstructed using a kinematical fit 
technique described in Ref . [55] . The incident neutron mo- 
mentum p n was calculated from the neutron velocity /3 n 
where the components p aiX and p niV were assumed to be 
zero. The set of kinematical variables x was chosen as 
x = {p n ',Pi, 6\, <j>i, ...,P3, #3, 4*3} where pt denotes the mo- 
mentum and 9i and (pi the polar and azimuth angle of the 
emitted particle i. 

Systematic errors in the kinematical fit due to system- 
atic shifts in the measured track parameters and veloci- 
ties were found to be much smaller than the experimen- 
tal resolution as indicated by the pull function pull{x) = 

^mcas _ x &ty A^T _ CT 2 mcaa; where ^meas &nd ^fit de _ 

note the measured and fitted values for the observables 
x, respectively. The resolutions for the pull functions were 
found to be between 0.9 and 1.0 with a good agreement be- 
tween experimental and Monte Carlo data indicating that 
the experimental errors had been reasonably estimated. 
The confidence level distribution CL shown in Fig. [6] for 
the data was found to be almost flat with a slight in- 
crease towards larger CL values. The strong enhancement 
at small CL values is caused by two kinds of events. About 
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Fig. 6. Confidence level distribution for data with full and 
empty target cell (cross hatched). 

half of these events originated from the target surround- 
ings (Fig. [6]). The other part are likely np — > pp7r~ events 
and background reactions of various kinds, e.g., np — > d7, 
np — ► d7r° and np — ► np7r° where electrons and positrons 
from the 7 conversion in the target surroundings could 
fake a fast pion. It was shown from Monte Carlo studies 
that about 8 % of np — > pp7r~ events give an enhance- 
ment at small CL values due to large multiple scattering, 
hadronic interactions and pion decay. A part of the np 
induced background reactions was already reduced by re- 
quiring a track intercept with the target cell and the cut on 
^vortex since the electron or positron track from the conver- 
sion process rarely matches with the proton or deuteron 
track. By rejecting all events with CL < 0.05 the recon- 
struction efficiency for the reaction np — > np7r° obtained 
from Monte Carlo was found to be 6 • 10~ 4 resulting in an 
estimated background contribution of 0.1 %. 

Although the cut at CL — 0.05 rejected a large amount of 
the background reactions in the LH2 target and the target 
surrondings, there were still background events from the 
Kapton walls of the target cell. As can be seen from Fig.lU 
these are 2-prong events where one track is associated with 
a pion and the other with a proton and are denoted 'pr- 
events in the following. They are mainly events from the 
quasi-free reactions on nuclei, nn — > d7r~ and in particular 
nn — > np7r~ , which have large cross sections due to the 
contribution of a±Q. A part of this background was rejected 
by the loose vertex cut whereas the main part could not 
be separated from the reaction np — > pp7r~ neither by the 
vertex reconstruction nor by the kinematical fit. 

For the reaction np — > pp7r~ , events of the 'p7r'-type oc- 
cured mainly when one proton stopped before reaching 
the DC8 drift chamber. The minimum momentum in the 
laboratory system for a final state proton from the reac- 
tion np — > pp7r _ decreases with increasing neutron energy. 
It reads 140 MeV/c at T n = 570 MeV. It was shown from 
Monte Carlo studies that protons with momenta less than 



190 MeV/c, when emitted from the central position in the 
target cell, stopped before reaching the hodoscope. At 570 
MeV, events of the 'p7r'-type from the reaction np — > pp7t _ 
contribute a few percent. This is confirmed by the data 
when comparing in Fig. 3] the 'p7r'-event distributions for 
data with full and empty target cell. Hence, it was decided 
to reject 'p7r'-events to reduce a large part of background 
while keeping most of the signal events. The still remain- 
ing background contribution from the target surroundings 
was of the order 4 % averaged over all neutron energies. It 
increased with decreasing neutron energy and was about 
8 % at 315 MeV. 




Fig. 7. The reconstructed vertex position z VC rtox for 2- and 3- 
prong events after all cuts for data with full and empty (grey) 
target cell. The background subtracted distribution is shown as 
full circles, Monte Carlo data as open circles. 

If the pion velocity and one proton velocity were almost 
of the same magnitude, the kinematical fit could find the 
wrong permutation. This happened in less than 5 % of 
the Monte Carlo generated events. However, due to the 
quite similar kinematical configurations for the correct 
and wrong reconstructed final state, this combinatorial 
background has a minor impact on the differential and 
integrated cross sections. 

The reconstructed z V ertex positions after all cuts are shown 
in Fig. [7] for the data with full and empty target cell. 
The background subtracted z VC rtex distribution is in good 
agreement with the expectation from the Monte Carlo 
simulation. 



4.7 Monte Carlo Simulations 

The experiment was simulated using the GEANT 3.21 
program package [56 incorporating all relevant materials 
of the experimental setup. The detector responses were 
implemented using detector efficiencies, resolutions and 
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calibrations as they were determined from the experimen- 
tal data, e.g., the effective light velocities in the scintilla- 
tors or the time-to-distance relations for the drift cham- 
bers. These responses were written in the same data stream 
format as the experimental data. 

Using 3.7 million phase-space distributed Monte-Carlo 
events, the reconstruction efficiency e pp7r ~ for the reaction 
np — > pp7r~ was determined as a function of the incoming 
neutron kinetic energy T n , the pion cm. angle 9* and the 
proton-proton invariant mass M pp : 



u pp7T 



PP 



(15) 



This function was a main input in the analysis in order 
to determine the invariant mass and pion angular distri- 
butions as well as the integrated cross sections. It was 
smoothed in order to minimize additional fluctuations in 
the data due to the limited Monte Carlo statistics. The 
average efficiency was of the order of 0.3 with a typical 
uncertainty of about 3%. The reconstruction efficiency 
e pp7r - shows a strong dependence as a function of M pp 
(see Sec l5.1[) . It drops towards the two-proton mass, since 
in this phase space region the two proton tracks are close 
together. In particular, for pions emitted in the backward 
direction, the value of e pp7r - at small proton-proton in- 
variant masses is only of the order of a few percent. The 
identical detector Monte-Carlo simulation was used to de- 
termine the reconstruction efficiency function e(T n , #*) for 
elastic np scattering. 



4.8 Elastic np scattering 

For detector calibration purposes and the neutron flux 
normalisation, np scattering events from the minimum 
bias trigger sample were selected. Assuming the elastic np 
scattering kinematics for events with one reconstructed 
track to be valid, the expected time-of-flight TOF^ xp for 
the incoming neutron was calculated using the measured 
track angle 8 P and the velocity f3f. Fig|5] shows the differ- 
ence ATOF n = TOF™P-TOF n + i- 19.75 ns, i = 0, 1, 
between the expected and measured neutron time-of-flight 
for data with full and empty target cell. Elastic np scatter- 
ing events kinematics gave values well-located at i ■ 19.75 
ns and were selected by the indicated cuts. For i — 0,1, the 
neutron time-of-flight corresponds to a neutron kinetic en- 
ergy above the pion production threshold. The main part 
of entries between the signal peaks is originating from re- 
actions on the Kapton walls and drift chamber materials. 
The remaining background part is likely due to inelastic 
reactions in the liquid hydrogen target like np — > d7r° or 
np — ► np7r° and was estimated to be about 1%. 



5 Results and discussion 

Differential cross sections as a function of the proton- 
proton invariant mass M pp and the pion cm. angle 9* and 
integrated cross sections were determined from the data. 
The M pp and angular distributions were subdivided in 
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Fig. 8. ATOF n distribution for data with full and empty 
(dashed line) target cell. Events within the indicated cuts (in 
grey) were accepted as elastic np scattering events. 



nine neutron energy bins where the first bin was between 
threshold and 330 MeV while the other bins were of 30 
MeV width. For the determination of the integrated cross 
sections, a finer binning of 10 MeV width was chosen in or- 
der to investigate the energy dependence of the cr np ^ pp7r - 
cross section in more detail. The differential cross sections 
da/dM P p and da/d9* have been normalised to yield the 
integrated cross section values at the corresponding neu- 
tron energies as given in Sec.| 



5.1 Invariant mass distributions 

The background subtracted and efficiency corrected M pp 
distributions for the nine neutron energy bins are shown 
in Fig.[9j The drawn errors include the statistical error as 
well as a 3% systematic uncertainty due to the efficiency 
function e pp7r - (T n , M pp ). The latter dominates for neutron 
energies above 345 MeV. For all energies the measured 
distributions differ significantly from the phase space ex- 
pectation. At higher energies, this is expected since the 
excitation of the A resonance contributes significantly to 
the production cross section an [35] . The strong final state 
interaction in the pp( 1 So) final state affects the Ss (en) 
and Sp (uoi) partial waves as already observed in the 
(Tn cross section pQ. Therefore, an enhancement at small 
Mpp values at least from the Ss partial wave is expected 
for neutron beam energies near threshold and is in fact 
observed in the data. For a qualitative understanding, 
the Mp P distribution from a Ss partial wave at T n — 
315 MeV is included in Fig. [9] where the detector reso- 
lution and the pp( 1 So) final state interaction have been 
taken into account. The final state interaction was calcu- 
lated in the effective range approximation [57j where the 
scattering length as = — 7.8098fm and the effective range 
r e ff = 2.767 fm were taken from Ref. 1581 . Distributions 




Fig. 9. Full dots: differential cross sections da/dM pp after background subtraction and efficiency correction. Histograms: nor- 
malised phase space expectations. Dotted lines: smoothed efficiency functions e -(T n ,M pp ) after integration overcosO^ where 
e pp7T - (T n , M pp ) = 1 corresponds to the maximum value of the ordinate. At 315 MeV, the M pp distribution from a Ss partial 
wave, scaled to the maximum of the measured data, is shown as a solid line. Data for the reaction np — > ppn~ at 443 MeV \4^j 
are shown as open circles. Dashed line: da/dM pp for the reaction pp — > pp7r° at 500 MeV \6 If '. 



produced by Sp partial waves from coi or Sd partial waves 
from (Tn would give a similar and even narrower signal. 

The ctoi contribution to the M pp distribution at T n = 
315 MeV has been estimated by a comparison to the re- 
action pp — > pp7T°. The Mpp distribution at 315 MeV, 
scaled by a factor of 2 due to ((2|), is shown in Fig. [10] to- 
gether with a measurement of the reaction pp — > pp7r° at 
T p = 310 MeV [50] which can be described by a combi- 
nation of Ss, Ps, Pp and Sd contributions [60] . Both dis- 
tributions differ in shape and magnitude which clearly in- 
dicates the presence of croi m the reaction np — > pp7r~ . 
The difference of 5 MeV in beam energies is slightly less 
than required for a comparison in the Q- or 77-scheme (« 7 
MeV). This has been partly compensated by scaling the 
M DD distribution at 315 MeV to the cross section value 



at 317 MeV by interpolating the cross section results for 
315 MeV and 325 MeV from Sec. 15. 31 The neutron energy 
bin width of about 40 MeV leads to entries for the reac- 
tion np — > pp7r~ beyond the kinematical limit for M pp 
from the reaction pp — > pp7r°- Both distributions were 
subtracted in order to extract the 001 contribution. Neg- 
ative values are likely due to the mismatch of the beam 
energies. By integrating the positive values, a ctqi cross 
section of about 0.9 /A> was found to be compared to a 
an value of about 4.3 ^b at T p = 310 MeV [60] . For a 
comparison, the M pp distribution from a Sp partial wave 
is shown. Again the detector resolution and the x So final 
state interaction, calculated in the effective range approx- 
imation, have been taken into account. The normalisation 
was chosen in order to yield the same ctqi cross section 
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value of 0.9 fib. The good agreement between the coi arL d 
the Sp distribution indicates that close to threshold the 
cross section aoi is mainly driven by Sp partial waves. 
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Fig. 10. Full boxes: da/dM pp at 315 MeV rescaled to give 
the integrated cross section at 317 MeV (see text) and multi- 
plied by a factor of two. Small dots: do jdM pp for the reaction 
pp — > pp7T° at 310 MeV [601. Open circles: difference between 
both distributions. Full dots: M pp distribution produced by a Sp 
partial wave taking into account detector resolution and final 
state interaction calculated in the effective range approxima- 
tion. 



Data of the experiment of Bachman et al. [ID] for the re- 
action np — >■ pp7r~ measured at 443 MeV are included in 
Fig. [9] Since the authors have published only relative cross 
sections, we normalised their distribution to our data. 
Their data also do not fit the phase space expectation but 
in addition show a strong deviation from our A/ pp distribu- 
tion at higher M pp values and as well near the two-proton 
mass. Due to the clear signal observed by Handler [35] at 
409 MeV, a 1 So enhancement at small M pp values has 
been expected |59j . The missing 1 So enhancement in the 
Bachman data was attributed to the poor invariant mass 
resolution which might dilute the signal [59 . However, it 
might be also due to an overestimation of the efficiency 
function in this particular phase space region. It should 
be noted that the result of Handler, given as a function 
of r = V%lv% max j i s m good agreement with the shape of 
our data [IS] . 

At 495 MeV, the M pp distribution is compared to the 
M pp distribution of the reaction pp — ► pp7r° at T p = 
500 MeV [ST] rescaled by a factor of 1 /2 (see Sect® . Since 
the Mpp distributions for both reactions differ significantly, 
one can conclude that croi gives also a large contribution 
to the reaction np — > pp7r~ at higher energies where the 
en cross section is already influenced by the excitation of 
the A resonance. 



5.2 Angular distributions 

The differential cross sections do^^^^- j dQ* for the 
nine neutron energies are shown in Fig.[TT] and Tab. [3] 
Corrections for reconstruction efficiency have been applied 
and background has been subtracted. For all energies, the 
angular distributions are anisotropic and at lower ener- 
gies, a pronounced f/b-asymmetry is observed. The solid 
curves are the results of a fit according to f [TTj) . The result- 
ing anisotropy parameters 6 np _> pp7r -, displayed in Fig. [12] 
show a strong dependence in T n with a maximum of 0.6 
at 375 MeV and a shallow minimum of 0.37 at 495 MeV. 
Compared to former experiments, the results show a sub- 
stantial improvement. Below 570 MeV, only two other ex- 
periments, at 409 MeV [2"5] and 435 MeVgDJ, have mea- 
sured over the full angular range. Our b np ^ pp7T - values are 
in quantitative contradiction to the data of Handler |28j , 
Kleinschmidt et al. |3T] and Bannwarth et al. [591 . 

The inclusive experiments of Kleinschmidt et al. [31] as 
well as Bannwarth et al. [39j did a pioneering work in es- 
tablishing the existence of coi but have some shortcom- 
ings. Both experiments were restricted in acceptance and 
their analyses relied on a model dependent extrapolation 
of the pion momentum spectra to small values. In the ex- 
periment of Kleinschmidt et al. [31J, positive charged pions 
were measured in the forward direction up to 8* ~ 30° 
only, which makes a reliable extraction of the angular 
distribution parameters difficult. In fact, the fits to the 
angular distributions were performed by setting the f/b- 
asymmctry parameter a identically zero. In the experi- 
ment of Bannwarth et al. [39] , positive and negative pions 
were measured at 6* = 90° and 6* = 166°. From these 
cross section values, the angular distribution parameters 
d\ and b np ^ pp7T - were determined using ([DJ). However, no 
explicit 7T^~ identification was performed and hence, sys- 
tematic errors might have been underestimated. 
The measured anisotropy parameters 6 np ^ ppw ~ are signif- 
icantly larger than those found in the reaction pp — > pp7r°, 
as shown in Fig. 1131 This is in qualitative agreement with 
most of the former findings, as can be seen from the com- 
parison of the results of older np experiments, shown in 
Fig-HH with the b pp ^ pp7r a values in FigfT3j This still holds 
when our data are compared only to the more recent 
6 P p^p p7 rO values of Rappenecker et al. [35] which are larger 
than the values of Dunaitsev et al. 162] and Stanislaus et 
al.JBTJ. 

The asymmetry parameters a shown in Fig. 1141 decreases 
monotonically with increasing energy which was already 
indicated by former experiments. However, quantitatively, 
the Bannwarth results [39J again deviate significantly from 
our data. Bachman et al. [ID] reported no significant f/b- 
asymmetry at 435 MeV but did not give a numerical 
value. We have fitted the angular distribution of Bach- 
man et al. [3DJ according to (TIT]) . We reproduced their 
value 6 np ^ ppw - = 0.47 ± 0.06 (Fig-UTS]) and found a = 
0.055 ±0.024 (Fig.[JJ|). 

The small value for a in the data of Bachman et al. [4Tj] 
might be attributed to the missing signal of the pp( 1 So) 




Fig. 11. Full dots: Differential cross sections dajdfl% after background subtraction and efficiency correction. Dotted lines: 
smoothed efficiency functions e pp7r - (T n , cos 6^) after integration over M pp where e pp7r - (T n , cos0£) = 1 corresponds to the max- 
imum value of the ordinate. Full lines: results of a fit according to 



final state interaction in their M pp spectrum. If their ef- 
ficiency function underestimated the contribution of this 
particular phase space region, the effect of the (I=0)-(I=1) 
interference between Ss and Sp partial waves would be 
suppressed. 

Both parameters, a and b ap ^ pp7r - , clearly indicate the 
significant contribution of the isoscalar production cross 
section ergi in the reaction np — > pp7T - , in particular in 
the energy range between 315 MeV and 400 MeV. 

The excitation of Pp(I=0) partial waves are expected to 
be suppressed at small energies due to their rf depen- 
dence and the anisotropies b pp ^, pp7! a were measured to 
be small for energies below 450 MeV. Since the partial 
wave 3 Si — s- 1 SoPi leads to a flat distribution in cos 8*, 
only the partial wave 3 Di — ^SoPi could explain the large 
&np->pp7i— values for energies below 450 MeV. The decreas- 
ing & n p->pp7i— values above about 375 MeV might be un- 



derstood in terms of increasing contributions from Ps par- 
tial waves in the 1=1 channel. The relative cross section 
from Ps partial waves in the reaction np — > pp7r _ was 
found to increase from 18% to 30% between T p = 325MeV 
and T p = 400MeV [64] . The b parameter tends to be con- 
stant above 450 MeV which can be understood in terms 
of increasing Pp partial waves from the isospin amplitude 
Mn being reflected in the rising of b pp ^ pp7r o observed by 
Rappenecker et al. [32] . 

The large f/b- asymmetries are likely due to the interfer- 
ence between the Sp partial waves from the isospin ampli- 
tude Mio and Ss from Mn. This is in agreement with the 
TPJUMF experiments [4l?ll4"4"] which have measured the Ss 
and Sp contributions in the reaction pn — > pp7r _ at very 
small Mpp values. With increasing energy, the interference 
signal vanishes rapidly. This can be understood by the rel- 
ative increase of the Ps and Pp partial waves in the 1=1 
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Table 3. Differential cross sections do/dfl* in ^ib/sr. 
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channel. In addition, due to the expected ry-dependence, 
the partial waves ratio Ss/Sp is supposed to drop with 
increasing beam energy. 



5.3 Integrated cross sections 



The Nnp^pp^- yields were obtained for each neutron en- 
ergy bin T n by integrating the events over M pp and cos 9* 
weighting each event by e"^- (T n , M pp , cos 9* ) , the inverse 
of the appropriate efficiency function value. The cross sec- 
tions were calculated from 



TV - 1 

-^np — >pp7r J- 



(16) 



where / is the lifetime factor of the data acquisition sys- 
tem. The time-integrated luminosity L is given by the 
product L = N n ■ F p , where F p denotes the number of 
protons per unit area and N n the time-integrated number 
of neutrons incident on the liquid hydrogen target. It was 
obtained by measuring the number of elastic np scatter- 
ing events in the minimum bias trigger sample at a given 
neutron energy and neutron cm. angle 9*: 



L(T n ,9* n ) 



1 



An* 



(in' 



f ■ P ■ e(T a ,e*) 



(17) 



The efficiency function e(T n , #*) for elastic np scattering 
was determined from the Monte Carlo simulation in steps 
of 10 MeV for T n and 10° for 6*, and P is the prescal- 
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Fig. 14. f/b- asymmetries compared to other np — > ppn~ ex- 
periments [28, 30, 38, 39,401. For the value quoted for the exper- 
iment of Bachman et al. \40\: see text. 



ing factor for the minimum bias trigger. The differen- 
tial cross sections dcr " rf p fj '" p were obtained using the par- 
tial wave analysis program SAID of Arndt et al. [65] . Fi- 
nally, for each neutron energy bin, the luminosity was av- 
eraged between 9* = 115° and 155°, the acceptance range 
where the efficiency function was approximately constant. 
Background contributions from the target surroundings 
for both, the A^p_, pp7r - (SecEH) and the N np ^ np (Sec Oj) 
yields, were subtracted using data taken with an empty 
target cell. In addition, the A np ^ np yields were corrected 
for the estimated background of about 1 % due to reac- 
tions in the LH 2 target (Sec. l4.8p . 

Compared to systematic errors statistical errors for the 
a np— >pp7i— cross sections were negligible for energies above 
325 MeV. A total systematic error of the order of 5.5% 
was estimated from three sources and these errors were 
added in quadrature: 

1 . the error of the efficiency function e pp7r - (T n , M pp , cos #* ) 
(3-4%), 

2. uncertainties in the experimental elastic np scattering 
yields and the efficiency function e(T n ,6!*) (3.5%), 

3. the uncertainty of the elastic np scattering differential 
cross section values (2.5%). 

The latter was estimated from the elastic cross sec- 
tion errors provided by the interactive partial wave 
analysis program SAID [65 (1%) and the difference of 
the cross section values between the Arndt analysis 
on one hand and the values given by the Nijmegen 
group [66] and Bystricky et al. [67] on the other hand. 
The Nijmegen analysis provides only values below 350 
MeV beam energy. The maximum observed deviations 
between the Nijmegen and the Arndt analysis are of 
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Table 4. Integrated cross sections o~ np ^ pp ^- in mb. 



the order 1%. The differential cross section values at 
9* = 135° given by the analysis of Bystricky et al. 
are significantly smaller than those given by the Arndt 
analysis. The deviation increases in magnitude from 
at 300 MeV up to -5% at 400 MeV and then decreases 
down to -1.5% at 560 MeV. 

It should be noted that neutron-proton elastic scat- 
tering cross section data taken at PSI [55] have been 
found to be significantly below the values provided by 
the Arndt analysis as well as the Bystricky analysis. 
The typical deviation is of the order —5% — 10%. 

The integrated cross section a np ^ pp7r - (Fig.[15l Tab. [4]) 
rises in the measured energy range by four orders of mag- 
nitude. The new data improve the knowledge of cr np ^p P 7r- 
for energies below 570 MeV substantially. 

The cross section measurements for the reaction np — * 
7r + X by Kleinschmidt et al. [31] which dominated the data 
between 480 MeV and 580 MeV so far are systematically 
below our data. This might be related to the above men- 
tioned acceptance cuts in their experiment. In addition, 
due to the different particle masses in the final states of 
the reactions np — * pp7r~ and np — * mi7r + , the cross sec- 
tion data can not be compared at the same beam energies. 
The T n values of Kleinschmidt et al. have to be lowered 
by about 6 MeV if the comparison would be performed 
in the ^-scheme resulting in a reduction of the observed 
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Fig. 15. Measured integrated cross section o~ n 
function of beam energy compared to other np — > 1N1N7T^ exper- 
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the Jiilich model \ia\70j . 
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Table 5. Fit parameters for <xn. 



more and more as the energy increases with a maximum 
deviation of a factor 2.5 at the largest energies. 



5.4 Extraction of er i 

From the measured a np ^ pp „- cross sections, Coi was cai ~ 
culated using The cross section an as a function of 
r\ was determined by a fourth order polynomial fit 



ou(»7) = 



(18) 



4=0 



deviation. Moreover, their data have to be corrected if 
improved cross section measurements are taken into ac- 
count. For the normalisation of their yields, the authors 
used cross section values of the reaction np — > d7r°. Those 
were determined by Hurster et al. [69 using the relation 
"np^dTr — §< 7 pp-»d7r+- Since this relation is only exact 
under the assumption of isospin invariance, several correc- 
tions due to the different particle masses and the Coulomb 
interaction in the pp- and d7r + system have to be applied. 
In consideration of the actual precise pp — > dn + data and 
all relevant corrections, the cr np _ >d-?r values of Hurster et 
al. [69] are too low by about 10% at 580 MeV and 20% 
at 480 MeV [55] • This results in an equivalent underesti- 
mation of the Cnp^Tr+x cross sections. If all these effects 
are taken into account, our cross sections and the Klein- 
schmidt data are found to be compatible. Nevertheless, the 
differences between our anisotropy parameters and those 
found by Kleinschmidt et al. still remain since the deter- 
mination of the anisotropy parameter does not depend on 
the flux normalisation. 

The most precise measurement at low energies, the data 
point of Handler [28] at 409 MeV, is significantly above 
our cross section values. This measurement averaged over 
a broad neutron energy spectrum where the mean neu- 
tron energy was determined by a maximum likelihood fit. 
It should be noted that the fit result for the neutron en- 
ergy spectrum lies slightly below the measured neutron 
energy spectrum [28] . As a consequence, the average neu- 
tron energy quoted by Handler might be underestimated. 

Our <T n p^ ppw - cross sections were also compared to the 
predictions of the Jiilich model [I0ll72] given as a full line 
in Fig. [15] At small neutron energies, the model overes- 
timates slightly the data points. Above T n « 320 MeV 
(r) ~ 0.5), the model underestimates our cross sections 



to the cross section data a pp ^ pp7I o of Meyer et al. [T], 
Bondar et al. [24], Rappenecker et al. [32], Stanislaus et 
al. [61] and Dunaitsev et al. [62] between T p — 285 MeV 
and T p = 572 MeV. The results of the fit parameters were 
mainly determined by the data sets of Meyer et al. [1] and 
Bondar et al. [24] for r\ < 0.6, and Rappenecker et al. [32] 
for -q > 1.15. The data of Dunaitsev et al. [62] and Stanis- 
laus et al. [61j which are filling the intermediate r\ range 
were needed for the proper convergence of the fit. The 
reduced x 2 value of \l = 2.165 with v = 63 degrees of 
freedom indicates some inconsistencies in the <7pp_> pp7r o 
data sets. The resulting fit parameters are listed in Tab. [5] 
and the result is shown in Fig. 1161 

The subtraction procedure according to ([2]) was per- 
formed within the rj- and the Q-scheme. The extracted 
(Toi-valties are shown in Fig. [15] as a function of rj. The 
errors were calculated from the propagated a np ^ pp7T - un- 
certainties. The different results for aoi found by the two 
methods are a measure of the systematic uncertainty in- 
duced by the different model assumptions. 

The obtained coi -values are reasonably described by a 
function of the form ooi <x 7 7 4 J a $ demonstrated by the 
solid line in Fig. 1161 except for the near-threshold region. 
The rf dependence strongly supports that urn is on ly car ~ 
ried by partial waves of the Sp type in the energy range 
below 570 MeV and confirms the interpretation of the M pp 
spectrum at 315 MeV. This finding is in contrast to the 
^9.i±2.4 dependence found by Kleinschmidt et al. [31] . Due 
to the large exponent, the authors concluded that mainly 
Pp partial waves should contribute to the isoscalar cross 
section. However, this strong rj dependence was caused by 
the above mentioned incorrect normalisation of the cross 
sections. 

The observed deviation of ooi from the rf behaviour near 
threshold might be addressed to the pp( 1 So) final state in- 
teraction. The modification of the cross section an near 
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integrated cross sections over four orders of magnitude 
have been determined resulting in a substantial improve- 
ment of the data compared to former measurements. A 
consistent picture of the reaction np — > pp7r _ has been 
found where all results establish the significant contri- 
bution of the isoscalar cross section aoi to the reaction 
np — > pp7r~ over the whole energy range. The determi- 
nation of the cross section aoi using the data of the re- 
action pp — + PPtt shows that aoi is mainly carried by 
Sp partial waves. The Jiilich model is able to describe the 
near-threshold a np ^ pp7r - cross section. However, with in- 
creasing energy, the prediction underestimates the data 
more and more where the main discrepancy is due to the 
description of the cross section an. The new data might 
provide also an interesting testing ground for calculations 
in the framework of chiral perturbation theory [22] . 



Fig. 16. The extracted cross section ooi in the rj- (full boxes) 
and Q-scheme (open circles) as a function of n. Drawn as a 
full line is a function proportional to -q 4 . For comparison, the 
parametrized an cross section (see text) is shown as a dotted 
line. The dashed line is aoi as predicted by the Jiilich model 
(see text). 

threshold due to the pp( 1 So) final state interaction is al- 
ready known [T] . However, the modification of the cross 
section an was mainly observed at much smaller values 

In the range 0.5 < i] < 0.9, the aoi cross section is at 
least of the same order as <Jn or even larger. This is the 
same region where the anisotropy parameters b np ^, pp7T - 
and the f/b-asymmetries are very large which supports 
their interpretation given in Sec. 15. 21 In this context, it is 
interesting to note that meson production models predict 
a dynamical suppression of the Ss partial wave [5|I72| in 
the same energy region. At higher values of r?, the resonant 
production with a Z\N intermediate state contributes more 
and more 35j and hence, an increases much stronger than 
(7oi. Nevertheless, aoi is found to contribute still about 
30% to the total a np ^ pp7r - cross section. 

Using the cross sections a np ^ PP Tr- and a pp ^ pp7T a pre- 
dicted by the Jiilich model [I0|I72|. the cross section ooi 
was calculated in the ry-scheme. As compared to a np ^ pp7r - 
(see FigfT5|). the cross section aoi is much better repro- 
duced by the Jiilich model (dashed line in Fig. [16]), al- 
though the excitation function shows a flatter slope than a 
rj 4 function resulting in a 50% underestimation of the cross 
section aoi at the largest energies. The underestimation 
of the <7np— I-PP7T- cross section at higher energies is mainly 
caused by the problem to describe the an cross section, 
in particular the higher partial waves Ps and Pp 64J. 



6 Conclusion 

We have measured the reaction np — > pp7r~ for neutron 
energies from threshold up to 570 MeV. Differential and 
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